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The rapid collapse of a polymer, due to external forces or changes in solvent, yields a long-lived 
‘crumpled globule.’ The conjectured fractal structure shaped by hierarchical collapse dynamics 
has proved difficult to establish, even with large simulations. To unravel this puzzle, we study a 
coarse-grained model of in-falling spherical blobs that coalesce upon contact. Distances between 
pairs of monomers are assigned upon their initial coalescence, and do not ‘equilibrate’ subsequently. 
Surprisingly, the model reproduces quantitatively the dependence of distance on segment length, 
suggesting that the slow approach to scaling is related to the wide distribution of blob sizes. 


The rapid collapse of a polymer into a dense globule, is 
a long-standing problem [lll~ll2j| . Such a collapse may be 
triggered by changes in solvent quality, causing the poly¬ 
mer to reduce its solvent-exposed surface area by forming 
a dense globule. A polymer may also be condensed by 
active forces, as in the rearrangement of DNA by pro¬ 
teins in the cell nucleus 0. The rapid collapse does not 
allow sufficient time for formation of topological entan¬ 
glements which abound in an equilibrated compact glob¬ 
ule 0-Q, .13]. It is suggested 0 that during collapse 

segments of the polymer initially condense to sphere-like 
‘blobs,’ which coarsen upon contact to form larger blobs. 
At any given time during the process the state is then as¬ 
sumed to be characterized by a single length scale 00 - 
ig]; e.g., the typical size of the blobs or the width of the 
tube connecting the blobs (see, e.g. Fig. mw 0o a 
central assumption is that when two blobs coalesce they 
remain more or less segregated within the newly formed 
structure. This is due to the slow relaxation processes 
within blobs, and due to topological constraints which 
prevent polymer segments forming the blobs from freely 
mixing - unlike a melt of independent polymer segments 
with open ends 0, 15l-[2Qj. The final configuration is thus 
predicted to be a constant density, self-similar, hierarchi¬ 
cal structure, known as the ‘crumpled’ or ‘fractal glob¬ 
ule’ 000. The end-to-end distance r m of segments 
of length m in the resulting globule is predicted to scale 
as r m ^ m l / d in d space dimensions (throughout the pa¬ 
per d = 3). This is in contrast to the equilibrium state 
reached at much later times 0 , where r m ^ ra 1 / 2 for 
small m, saturating at the globule size r max = N 1 ^, 
where N is the length of the polymer 0- 


These predictions have been tested in several simula¬ 
tions of polymer compaction 0000, which gener¬ 
ally confirm that the rapidly collapsed state is not en¬ 
tangled, and indeed different from the equilibrium glob¬ 
ule. However, they do not agree upon its fractal na¬ 
ture. In particular, the expected scaling r m ^ m l ^ d has 
not been conclusively confirmed, even with the largest 
size simulations (recently extended to polymers of up to 
250,000 monomers |6j, Hlj). This implies either a very 
large crossover scale before fractal behavior to clearly 
manifested, or else that the collapsed state is not strictly 



FIG. 1. The Coalescence Model combines a simulation of co¬ 
alescing spherical drops (b,d,e), with estimates for distances 
between monomer pairs in final structure. MD simulations 
(a,c) with initial conditions identical to (b,d) are shown. Ini¬ 
tial conditions are SAW for (a,b), RW for (c,d), and ID for 
(e). All panels show parts of the full systems. Colors show 
position of monomers along the polymer, or average position 
for drops. 


self-similar. The large finite-size effects are at least partly 
explained by partial mixing of short polymer segments, 
for which the topological constraints do not apply 0|l5- 
I 20 L HI]. More generally, various protocols for con- 
structing crumpled globules have been suggested 00 
(i. ijj, [ 23 ], and it is not clear if the different procedures 
yield the same state. Settling these issues seems to re¬ 
quire even larger simulations, or new theoretical insights. 


To this end, we propose a coarse-grained model for 
the crumpled globule formed by polymer collapse, fo¬ 
cusing here on active compression. The evolution of 
blobs is modeled by drop coalescence, while distances 
between monomers in the final crumpled state are as¬ 
signed without keeping track of their individual posi¬ 
tions. This builds upon the topological segregation of 
blobs, and the slowness of subsequent internal rearrange- 
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merits, and should apply at scales beyond which these 
conditions hold. The results highlight a different, dy¬ 
namical source of slow convergence to a self-similar state, 
and draw connections to the physics of coagulation and 
drop-coalescence. In particular, the assumption of single 
length-scale during the collapse may need to be amended, 
at least when the collapse is due to active compression. 

In spirit of the blob theory 0 , the basic entity in our 
model is a ‘drop,’ an abstraction of the blob. A drop 
is a uniform density sphere which contains a subset of 
the monomers from the original polymer, without ex¬ 
plicit position assignments within the drop. Initially, ev¬ 
ery monomer is a drop. Drops move (as detailed below) 
and coalesce into larger drops immediately upon colli¬ 
sion, i.e. as soon as the spheres overlap in space. Drops 
do not break into smaller drops. Coalescence conserves 
volume; drops of volume V\ J V 2 forming a drop of vol¬ 
ume v\ + V 2 - (With the monomer volume set to unity, 
the drop volume equals the number of monomers it con¬ 
tains.) The new drop is centered at the center of mass 
of the coalesced drops. The process terminates when all 
drops have merged to a single sphere. Example coales¬ 
cence runs are compared with corresponding molecular 
dynamics (MD) runs in Fig. [lj 

Our main interest is the structure of the final collapsed 
state, as captured by the distances {r m }. However, in the 
coarse-grained drop coalescence model we do not keep 
track of the internal structure of a drop. Instead we as¬ 
sign distances-estimates to all pairs of monomers within 
a drop. Guided by the slow internal rearrangements in 
the blob picture, a distance-estimate is assigned when 
a pair of monomers first comes together in a coalescence 
event, and is not changed thereafter. Upon coalescence of 
drops of volumes v\ and V 2 , v\ XV 2 pairs of monomers, one 
from each drop, are assigned a distance-estimate. This 
distance is of order (v\ +V 2 ) 1 ^ d , the linear size of the 
new drop, which can be sampled from any distribution. 
As the results are highly insensitive to this choice, we 
simply assign all V\V 2 pairs the ‘distance’ (v\ + t^) 1 ^. 
The distance-estimates satisfy triangle inequalities, and 
importantly are ultrametric, the tree structure reflecting 
the hierarchy of the collapse process [Hj. At the end, 
when all monomers belong to a single drop, all monomer 
pairs have been assigned a distance-estimate. 

To fully define the model, it remains to prescribe the 
drops’ motion between coalescence events, as well as the 
initial conditions. Here we focus on a simple dynamics for 
active compression in which the drop velocity is propor¬ 
tional to its distance from the origin, i.e. dR a /dt = — R a , 
where R a is the position of the drop center. This cor¬ 
responds to over-damped motion in a harmonic poten¬ 
tial. It can also be viewed as a uniform compression of 
space, as the distance R a p = R a — Rp between two drops 
that have not yet coalesced changes as dR a p/dt = —R a f3, 
so all relative distances shrink by the same factor per 


unit time. Apart from coalescence, there is no addi¬ 
tional interaction between drops. In particular, we do 
not impose polymeric bond interactions between sequen¬ 
tial monomers. 

The initial conditions can be any configuration of the 
monomer positions. Here we use random-walk (RW), 
self-avoiding walk (SAW), and ‘one-dimensional’ (ID) 
initial conditions. For the latter, monomers are posi¬ 
tioned along a line as Rj = jx + rjj, where x is the unit 
vector in the x-direction, and r]j are independent Gaus¬ 
sian variables in d dimensions with (rjj) = 0, (r/J) = c 2 , 
c a number of order of drop radius. Open polymers are 
used in all initial conditions. The SAW initial-conditions 
are arguably the most natural for a polymer in a good 
solvent, but other initial conditions help in obtaining ad¬ 
ditional insight. In particular, ID initial conditions are 
used only in the coalescence process, as the strong unidi¬ 
rectional compression will cause large internal rearrange¬ 
ments after blobs have formed within MD. 

The coalescence model is compared with MD simula¬ 
tions where the polymer is composed of monomers with 
standard polymeric bond attraction forces F nn , and re¬ 
pulsion F rep between all monomer pairs, see Supplemen¬ 
tal Material. Together with the same external force 
as in the coalescence model, the monomers evolve as 
dRi/dt = —Ri + F nn + F rep . The MD runs are termi¬ 
nated when the polymer size stops decreasing. In both 
the coalescence and MD simulations we do not add noise 
to the dynamics, as the entire collapse process is fast (of 
order In TV, see below). Noise may further increase re¬ 
arrangement processes in the MD simulations, which we 
try to minimize here. A long-lived, metastable state is 
obtained, due to the time-scale separation between the 
fast collapse process and further rearrangements which 
are much slower. (As noted above, other ‘crumpled glob¬ 
ule’ metastable states may also exist jD-EI 22], at later 
times.) 

Now, let dif al be the distance-estimate between 
monomers in the coalescence model, and d^ D = 
\\Ri — Rj\\ the corresponding distance between the 
monomers in the final state of the MD run. In what 
follows we apply the same analysis to d^ al , df!j D , ref- 
ered jointly as dij. We define the normalized end-to- 


end distance r m = Cy (d? J -)|*_^| s=m , where the average 

( m ")\i-j\=m over m o n o m er pairs separated along 
the polymer by m monomers, as well as over repeated 
runs of the two models, with initial conditions resam¬ 
pled for each run. The normalization C is chosen so that 
E™ jv(N-T)/ 2 r m ~ 1 - ( N - m is simply the number 
of pairs that are m apart in an open polymer of length 
TV). In this way the MD and coalescence models can be 
compared without any fitting parameters. 

Figure [2fa) compares r m from the coalescence model 
with MD simulations for SAW initial conditions. MD 
results are shown for different polymer lengths TV. Strik- 
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FIG. 2. Comparison of Coalescence model and MD results, 
(a) Normalized end-to-end distance r m , with SAW initial 
conditions. Inset: r m / ( m/N ) 1 ^ 3 for the same V m • (b) 
rm/ (m/N) 1 / 3 for RW initial conditions, (c) r m / (m/N) 1 / 3 
for Coalescence model with different initial conditions, and 
system sizes, N. 


ingly, the MD and coalescence models agree quantita¬ 
tively over the entire range where the MD results of 
different polymer sizes collapse, for 10 -3 < m/N < 1. 
In the small m regime, where MD simulations with dif¬ 
ferent N do not collapse, they show a different behav¬ 
ior (consistent with a random walk) not present in the 
coalescence model. This demonstrates that the coales¬ 
cence process indeed represents a coarse-grained model, 
capturing large-scale behavior while incorporating differ¬ 
ent microscopic details. The absence of polymeric bonds 
in the coalescence model removes short-scale constraints 
that (as supported by the MD simulations) do not effect 
large-scale properties through the distance assignment. 

The m 1 / 3 trend-line in Fig. [2fa) indicates that the ex¬ 
pected scaling is not present at the tested system sizes 
N < 25 x 10 3 . To more carefully assess self-similarity, we 
study finite-size scaling in the standard form 

r m = (m/N) 1/3 f (m/N) , (1) 

expected to hold for 1 m. (The additional TV -1 / 3 
reflects the choice of normalization.) The scaling func¬ 
tion / (x) should go to a constant for x 0, such that 
r m oc (m/N) 1 / 3 for 1 m N. Data for both MD and 
coalescence models, for different initial conditions and 
system sizes, are summarized in Fig. [2j The difference 
between MD and coalescence results for the RW initial 
conditions is larger than for the SAW. Importantly, for 
both SAW and RW initial conditions, and for both mod¬ 
els, the expected condition, f (m/N 0) = const , is not 
seen clearly for tested values of TV; in all cases the maxi- 
mum of f (m/N) = r m / (m/N) 7 appears to grow with 
increasing N. This effect is almost absent for the ID ini¬ 
tial conditions, and largest in the RW case. Since the coa¬ 
lescence model includes only a minimal set of microscopic 
details, it is surprising to see such a slow approach to the 


expected scaling behavior. Understanding this trend in 
the coalescence model should provide insights into the 
more complex case of the collapsing polymer. 

While simplified, a full understanding of the coales¬ 
cence model - including the distribution of distance- 
estimates assigned as a function of the time t and sep¬ 
aration m - is still difficult. Fortunately, some insights 
regarding scaling (or lack thereof) can be gained by ex¬ 
amining the distribution of drop volumes as a function 
of time, even without making reference to the distance- 
estimate assignments. In particular, the volume distri¬ 
butions already show a slow approach to scaling. 

Let p t ( v ) be the distribution of drop sizes at time t. 
Volume conservation implies that f dvvp t (v) = TV, while 
f dvpt (v) is the number of drops at time t (averaged over 
repeated runs). Within the scaling picture, the dynamics 
of the collapse is described by a single typical drop size 
as a function of time, v(t), when 1 <C v (t) <C N. For 
example, in the tube picture v 1 / d may be the thickness of 
the tube. In such a scaling regime, p t (v) should depend 
on time only through v (t), as 


pt («) 


N 


V) 


( 2 ) 


The factor v 2 ensures that / dvvpt (v) remains con- 
stant. We measure the typical size via v (t) = 

(E„ v n)/E„ v n = (E n V l)l N i where the sums run 
over all drops at time £, and the average (• • •) is over 
initial conditions. This definition, standard in drop co¬ 
alescence/coagulation literature 23-[27j. addresses possi¬ 
ble divergences of p t (v) at small volumes (see below). To 
test Eq. [2j we plot </> t (v) = N _1 v ( t) 2 p t (v), against the 
normalized volume v/v (t), and check for data collapse at 
different values of t. 

The distributions, depicted in Figs. [3ja,b) for ID and 
RW initial conditions, respectively, are quite different. 
The one in (a) (ID initial conditions) is concentrated in 
the region 0.4 < v/v < 1.6, and strongly suppressed out¬ 
side this interval. Thus, at any given time all volumes 
are of the same order, with a ratio of about 4 between the 
largest and smallest. The distributions at different times 
collapse nicely when plotted against v/v (t). In contrast, 
the distributions for RW initial conditions in Fig. (U^b) 
are very wide (note the log-log scale), with possibly a 
diverging tail at small volumes, v/v 1. (The dashed 
line, x _1 , is included as an indication of such potential 
divergence.) Moreover, the distributions fail to collapse 
in this tail. With SAW initial conditions, Figs. (3fc), the 
results appear to be intermediate between the above two 
cases, with a distribution that is finite at v = 0 (at least 
for tested N). The simulations in Figs.[Ua,b,c) were car¬ 
ried out with N = 2.5xl0 4 , 5x10 4 , and 10 5 respectively, 
to allow for better scaling in Figs. Etb,c). 

The evolution of v (t) is depicted in Fig. [3jd). For 
ID and SAW initial conditions v (t) oc e wt is a good fit 
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FIG. 3. Drop volume distributions in the Coalescence model, 
(a) Id initial conditions. Dashed line: sequential collisions 
theory, (b) RW initial conditions, dashed line: x~ x trendline, 
(c) SAW initial conditions, (d) v(t) for the three models. 
Dashed lines: e wt for different cases (see text). 


at intermediate times, where w = du ^_ 1 , with z/q de¬ 
scribing the scaling ~ m u ° in the initial condition. 
This form is explained as follows: At time £, segments 
of initial length v (t) form drops of volume v (t ), with 
diameter [i v (t)] 1 ^ d . The initial end-to-end distance of 
these segments scales as [v . If this reduction in 
length follows the general exponential approach of free 
monomers [28|, then \v (t)] u ° e~ l — [^(t)] 1 ^, leading to 
the above relation [29|. For RW initial conditions the 
growth of v (t) follows this form for a narrower interval, 
with wide crossover regions, probably related to the lack 
of scaling observed in Fig. [3^b). 

The difference between the three volume distributions 
in Fig. 02 can be directly observed in the coalescence and 
MD runs in Fig. [TJ note especially the large number of 
small drops for RW initial conditions in Fig. QJd). In 
the MD starting from the same initial conditions, these 
manifest as open segments of the polymer, alongside large 
condensed regions, see Fig. [ljc); the thickness of a puta¬ 
tive ‘tube’ is highly uneven. We now present theoretical 
approaches to analogous problems leading to the two very 
different volume distributions obtained in Figs. (3ja,b). 


Broad distributions, with power-law tails at small vol¬ 
umes, as in Fig. [3jb) for RW initial conditions, appear 
in related problems such as diffusion-limited aggrega¬ 
tion [3jo|, and drop coalescence 24, [25]. Heterogeneous 
drop coalescence, where drops grow but no new drops 
are added, is perhaps most similar to our model. There, 
the distribution of drop sizes is highly poly-disperse when 
v (t) (in our notation) grows exponentially in time [24|. 
Aggregating systems are commonly studied through the 


(mean-field) Smoluchowski equation 0113: 

d t p(v) = \ J dv'Jpy , v -v' ~ J dv'Jt-v 


(3) 


The first term is the change in density p t (v) due to 
creation of drops of size v, while the second term de¬ 
scribes their removal due to coalescence. In the current 
context, this approximation postulates a rate Jt-vi,v 2 = 
K Vl :V2 Pt ( vi ) pt (V 2 ) for collisions between drops of vol¬ 
umes Vi,V 2 , with a kernel K Vl , V2 depending only on the 
volumes of the coalescing drops. Here, rather than ex¬ 
plicitly constructing an approximate kernel, we refer to 
extensive ly s tudied scaling solutions of the Smoluchowski 
equation [26|, (27j. In particular, for homogeneous kernels 
such that K aVl ,av 2 = a x K VuV2 for a > 0, it is known that 
when v (t) oc e wt , as in our case, the scaling function 
4>{x) in Eq. [2] has the following properties: It is strongly 
(exponentially) decaying for x 1, and has a diverging 
power-law tail (x ) ^ x~ T for 1 x, with 1 < r < 2 
(the value of r depends on the kernel). For example, the 
kernel K Vl ^ V2 = v\ + admits an exact scaling solution 
<fi (x) = -==x~ 3 / 2 e- X//2 with v = e 2t , so that r = 3/2. 

Interestingly, a slow approach to the asymptotic scal¬ 
ing is well documented for certain cases with such small- 
volume tails [13, MM • The nature of this slow approach 
is still not well-understood, and might be sensitive to the 
kernel form and initial conditions. It is intriguing to spec¬ 
ulate on its relation to the present problem. 

In the case of ID initial conditions in Fig. [3fa), essen¬ 
tially all drops at a given time have volumes of the order 
of v(t). Unlike RW initial conditions, all collisions here 
are sequential, i.e. between drops comprised of adjacent 
segments along the polymer. Here geometry matters: 
after a collision gaps are formed on both sides of the 
newly constructed drop, greatly reducing its probability 
of coalescing again before other drops have time to grow. 
Smaller drops leave smaller gaps, and have an increased 
probability of additional collisions. These effects are dis¬ 
cussed in a quantitative way in the Supplemental Mate¬ 
rial. An approximate evolution equation is derived that 
admits the scaling solution shown in Fig. (3ja) (dashed 
line), which strongly decays outside a narrow interval of 
v/v, just like the results from the full coalescence model. 

The coalescence model proposed here does away with 
several microscopic details present in the collapsing poly¬ 
mer that could delay the asymptotic approach to scaling. 
The lack of simple scaling in coalescence thus points to 
deeper problems with the simple model of hierarchical 
collapse, in particular in the assumption of a sharply 
defined blob-scale. The coalescence model is interest¬ 
ing in its own right; it is closely related to widely stud¬ 
ied problems of coarsening of growing water drops, but 
differs in the initial (polymeric) distribution of droplets. 
As demonstrated, different such initial conditions lead to 
widely dissimilar probability distributions. Probing the 
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role of fractal initial conditions in coalescence problems 
should thus be worthy of further exploration. 
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SUPPLEMENTAL MATERIAL 

MOLECULAR DYNAMICS (MD) SIMULATION 
DETAILS 


In the MD simulations, the polymer is composed of 
monomers with standard polymeric bond forces F nn , and 
repulsion F rep between all monomer pairs. F nn , F rep are 
given by the gradients F nn = -W nn , F rep = -W rep , 
where Ilf 


U nn (r) = | ^ 

Urep (r) = \ 4£ 


Ik/ 2 In 


1 - ( r/lmax ) 


00 

vl2 / / \6 


r < ^ 
r > L 


(a/r) - (a/r) J + C (r cu t) r < r cu t 

0 r > r cu t 


with k = 30, l max = 1.5, e = a = 1 (H, and r cut = 2 1 / 6 
(which makes U rep purely repulsive). C (r cut ) is chosen 
such that U rep (r = r cut ) = 0. 

Together with the same external force as in the coales¬ 
cence model, the monomers evolve as dRi/dt = — Ri + 
F nn + F rep . We use first-order integration in time. The 
MD runs are terminated when the polymer size stops 
decreasing. Of course, there may be further internal re¬ 
arrangements due to reptation, but on much longer time- 
scales. The final states are nearly spherical, with at most 
a few percent ellipticity. 


THEORY FOR DROP EVOLUTION ON A LINE 


taken to be one). If these segments form drops, they 
will meet when (^ d + e 1 ~ ( - Vl 1 2 v2 ± cyjv\ + v ^). 

The resulting time frame is narrow when 1. Small 

drops are rapidly swept away as they coalesce with typi¬ 
cal sized drops at earlier times than larger drops. 

The above argument can be turned into an evolution 
equation for p+ (v). For other approaches to a related 
problem, see [2, |3[. We postulate a collision rate 

Jt\V 1,V2 = K-Vi ,V2',tPt (M ) Pt (^ 2 ) 5 (Si) 

as in Equation 3 in the Main Text, but with a kernel 
K Vl ,v 2 ]t that explicitly depends on time (therefore this 
is not a Smoluchowski equation). The kernel represents 
the probability per unit time for collision of drops t’l, ^ 2 , 
given that they have not yet coalesced. It corresponds 
to (■ dS 2 /dt ) /S 2 , where ^(tj^i,^) is the survival prob¬ 
ability of the two drops till time t. We approximate S 2 
by the probability that two drops formed by a random 
pair of consecutive segments do not overlap at time £, 
i.e., Sxoe~ f > v\^ d + 


1 V 2 

S 2 (t-,vi,v 2 ) = / e~^dr/ (S2) 

\/27rC JG(v-\_,V2,t) 

where 

e* (vl /d +vl /d ) - v -^ 

G(vi,v 2 ,t) = --- , (S3) 

so that the kernel takes the form 


Here we further discuss the evolution of drop sizes for 
ID initial conditions, and describe a theory that accounts 
for the distribution of volumes in Fig. 3(a). 

A discussed in the Main Text, for ID initial condi¬ 
tions essentially all drops at a given time have volumes 
of the order of v(t), see Fig. 3(a). All collisions here 
are sequential, i.e. between drops comprised of adjacent 
segments along the polymer. Here geometry matters: 
after a collision gaps are formed on both sides of the 
newly constructed drop, greatly reducing its probability 
of coalescing again before other drops have time to grow. 
Smaller drops leave smaller gaps, and have an increased 
probability of additional collisions. 

To quantify this effect, note that a drop moves accord¬ 
ing to R(t) = where Ro is the center of mass of 

the monomers’ initial positions. Consider two sequential 
segments, containing v\ and V 2 monomers respectively. If 
the segments are chosen at random (while still being adja¬ 
cent), the distance between their centers is Gaussian 
distributed, with average Vl + V2 and variance c(v\ +^ 2 ) 
(the average distance between neighboring monomers is 


/^i 2 (—1 2 

_ e~~ d t G _ e~~ d t G 
erfcG 


(S4) 


Equation 3 with this kernel indeed has a scaling solution 
with v oc e wt , as can be verified by direct substitution. 
Numerically solving Eq. 3 with this kernel, one obtains 
a scaling collapse with v quickly approaching e wt , and 
a scaling form 0 (|) shown in Fig. 3(a) (dashed line), 
which strongly decays outside a narrow interval of v/v, 
just like the results from the full coalescence model. 
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